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Abstract 

We review the q-deformed spin network approach to Topological Quantum Field 
Theory and apply these methods to produce unitary representations of the braid 
groups that are dense in the unitary groups. Our methods are rooted in the 
bracket state sum model for the Jones polynomial. We give our results for a 
large class of representations based on values for the bracket polynomial that 
are roots of unity. We make a separate and self-contained study of the quan- 
tum universal Fibonacci model in this framework. We apply our results to 
give quantum algorithms for the computation of the colored Jones polynomi- 
als for knots and links, and the Witten-Reshetikhin- Turaev invariant of three 
manifolds. 



Introduction 



This paper describes the background for topological quantum computing in 

terms of Temperley-Lieb recoupling theory. This is a recoupling theory that 
generalizes standard angular momentum recoupling theory, generalizes the 
Penrose theory of spin networks and is inherently topological. Temperley-Lieb 
recoupling Theory is based on the bracket polynomial model [34, 41] for the 
Jones polynomial. It is built in terms of diagrammatic combinatorial topology. 
The same structure can be explained in terms of the SU{2)q quantum group, 
and has relationships with functional integration and Witten's approach to 
topological quantum field theory. Nevertheless, the approach given here will 
be unrelentingly elementary. Elementary, does not necessarily mean simple. 
In this case an architecture is built from simple beginnings and this archic- 
tecture and its recoupling language can be applied to many things including, 
e.g. colored Jones polynomials, Witten-Reshetikhin-Turaev invariants of three 
manifolds, topological quantum field theory and quantum computing. 

In quantum computing, the application is most interesting because the 
simplest non-trivial example of this Temperley-Lieb recoupling Theory gives 
the so-called Fibonacci model. The recoupling theory yields representations 
of the Artin Braid group into unitary groups U{n) where n is a Fibonacci 
number. These representations are dense in the unitary group, and can be 
used to model quantum computation universally in terms of representations of 
the braid group. Hence the term: topological quantum computation. 

In this paper, we outline the basics of the TL- Recoupling Theory, and 
show explicitly how the Fibonacci model arises from it. The diagrammatic 
computations in the section 9 are completely self-contained and can be used 

by a reader who has just learned the bracket polynomial, and wants to see how 
these dense unitary braid group representations arise from it. The outline of 
the parts of this paper is give below. 

1. Knots and Braids 

2. Quantum Mechanics and Quantum Computation 

3. SU{2) Representations of the Artin Braid Group 

4. The Bracket Polynomial and the Jones Polynomial 

5. Quantum Topology, Cobordism Categories, Temperley-Lieb Algebra and 
Topological Quantum Field Theory 
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6. Braiding and Topological Quantum Field Theory 

7. Spin Networks and Temperley-Lieb Recoupling Theory 

8. Fibonacci Particles 

9. The Fibonacci Recoupling Model 

10. Quantum Computation of Colored Jones Polynomials and the Witten- 
Reshetikhin-Turaev Invariant 

We should point out that most of the results in this paper are either new, 
or are new points of view on known results. The material on SU (2) represen- 
tations of the Artin braid group is new, and the relationship of this material 
to the recoupling theory is new. See Theorem 1 in Section 3. The treatment 
of elementary cobordism categories is well-known, but new in the context of 
quantum information theory. The reformulation of Temperley-Lieb recoupling 
theory for the purpose of producing unitary braid group representations is 
new for quantum information theory, and directly related to much of the re- 
cent work of Freedman and his collaborators. In Theorem 2 in Section 7 we 
give a general method to obtain unitary representations at roots of unity us- 
ing the recoupling theory. The treatment of the Fibonacci model in terms 
of two-strand recoupling theory is new and at the same time, the most ele- 
mentary non-trivial example of the recoupling theory. The models in section 
10 for quantum computation of colored Jones polynomials and for quantum 
computation of the Witten-Reshetikhin-Turaev invariant are new. They take 
a particularly simple aspect in this context. 

Here is a very condensed presentation of how unitary representations of the 
braid group are constructed via topological quantum field theoretic methods. 
One has a mathematical particle with label P that can interact with itself to 
produce either itself labeled P or itself with the null label *. When * interacts 
with P the result is always P. When * interacts with * the result is always *. 
One considers process spaces where a row of particles labeled P can successively 
interact subject to the restriction that the end result is P. For example the 
space K[(a6)c] denotes the space of interactions of three particles labeled P. 
The particles are placed in the positions a, 6, c. Thus we begin with {PP)P. 
In a typical sequence of interactions, the first two P's interact to produce a *, 
and the * interacts with P to produce P. 

{PP)P — > (*)P — > P. 
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In another possibility, the first two P's interact to produce a P, and the P 
interacts with P to produce P. 

{PP)P — >{P)P — > P. 

It follows from this analysis that the space of linear combinations of processes 
K[(a6)c] is two dimensional. The two processes we have just described can 
be taken to be the qubit basis for this space. One obtains a representation 
of the three strand Artin braid group on V[(a6)c] by assigning appropriate 
phase changes to each of the generating processes. One can think of these 
phases as corresponding to the interchange of the particles labeled a and b in 
the association [ab)c. The other operator for this representation corresponds 
to the interchange of b and c. This interchange is accomplished by a unitary 
change of basis mapping 

F : V[{ab)c] — > V[a{bc)]. 

If 

A : V[{ab)c\ — > V[{ba)c\ 

is the first braiding operator (corresponding to an interchange of the first two 
particles in the association) then the second operator 

B : V[{ab)c] — > V[{ac)b] 

is accomplished via the formula B = F^^RF where the R in this formula acts 
in the second vector space V[a(6c)] to apply the phases for the interchange of 
b and c. These issues are illustrated in Figure 1, where the parenthesization 
of the particles is indicated by circles and by also by trees. The trees can be 
taken to indicate patterns of particle interaction, where two particles interact 
at the branch of a binary tree to produce the particle product at the root. See 
also Figure 28 for an illustration of the braiding B = F~^RF 

In this scheme, vector spaces corresponding to associated strings of particle 
interactions are interrelated by recoupling transformations that generalize the 
mapping F indicated above. A full representation of the Artin braid group 
on each space is defined in terms of the local interchange phase gates and the 
recoupling transformations. These gates and transformations have to satisfy 
a number of identities in order to produce a well-defined representation of the 
braid group. These identities were discovered originally in relation to topolog- 
ical quantum field theory. In our approach the structure of phase gates and 
recoupling transformations arise naturally from the structure of the bracket 
model for the Jones polynomial. Thus we obtain a knot-theoretic basis for 
topological quantum computing. 
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Figure 1 - Braiding Anyons. 

Aspects of the quantum Hall effect are related to topological quantum 
field theory [77, 23, 14, 15], where, in two dimensional space, the braiding of 
quasi-particlcs or collective excitations leads to non-trival representations of 
the Artin braid group. Such particles are called Anyons. It is hoped that the 
mathematics we explain here will form the bridge between theoretical models 
of anyons and their applications to quantum computing. 
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1 Knots and Braids 



The purpose of this section is to give a quick introduction to the diagrammatic 
theory of knots, links and braids. A knot is an embedding of a circle in three- 
dimensional space, taken up to ambient isotopy. The problem of deciding 
whether two knots are isotopic is an example of a placement problem, a problem 
of studying the topological forms that can be made by placing one space inside 
another. In the case of knot theory we consider the placements of a circle inside 
three dimensional space. There are many applications of the theory of knots. 
Topology is a background for the physical structure of real knots made from 
rope of cable. As a result, the field of practical knot tying is a field of applied 
topology that existed well before the mathematical discipline of topology arose. 
Then again long molecules such as rubber molecules and DNA molecules can 
be knotted and linked. There have been a number of intense applications of 
knot theory to the study of DNA [72] and to polymer physics [55] . Knot theory 
is closely related to theoretical physics as well with applications in quantum 
gravity [76, 71, 50] and many applications of ideas in physics to the topological 
structure of knots themselves [41]. 



Quantum topology is the study and invention of topological invariants via 
the use of analogies and techniques from mathematical physics. Many invari- 
ants such as the Jones polynomial are constructed via partition functions and 
generalized quantum amplitudes. As a result, one expects to see relationships 
between knot theory and physics. In this paper we will study how knot the- 
ory can be used to produce unitary representations of the braid group. Such 
representations can play a fundamental role in quantum computing. 
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Figure 2 - A knot diagram. 




Figure 3 - The Reidemeister Moves. 

That is, two knots are regarded as equivalent if one embedding can be obtained 
from the other through a continuous family of embeddings of circles in three- 
space. A link is an embedding of a disjoiint collection of circles, taken up to 
ambient isotopy. Figure 2 illustrates a diagram for a knot. The diagram is 
regarded both as a schematic picture of the knot, and as a plane graph with 
extra structure at the nodes (indicating how the curve of the knot passes over 
or under itself by standard pictorial conventions) . 
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Braid Generators 



S2 = S2 S2 



S-j Sg = SrjS 



3 ^1 



Figure 4 - Braid Generators. 



Ambient isotopy is mathematically the same as the equivalence relation 
generated on diagrams by the Reidemeister moves. These moves are illus- 
trated in Figure 3. Each move is performed on a local part of the diagram 
that is topologically identical to the part of the diagram illustrated in this 
figure (these figures are representative examples of the types of Reidemeister 
moves) without changing the rest of the diagram. The Reidemeister moves are 
useful in doing combinatorial topology with knots and links, notably in work- 
ing out the behaviour of knot invariants. A knot invariant is a function defined 
from knots and links to some other mathematical object (such as groups or 
polynomials or numbers) such that equivalent diagrams are mapped to equiv- 
alent objects (isomorphic groups, identical polynomials, identical numbers). 
The Reidemeister moves arc of great use for analyzing the structure of knot 
invariants and they are closely related to the Artin Braid Group, which we 
discuss below. 
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Figure 5 - Closing Braids to form knots and links. 




b CL(b) 

Figure 6 - Borromean Rings as a Braid Closure. 

Another significant structure related to knots and links is the Artin Braid 
Group. A braid is an embedding of a collection of strands that have their ends 
top and bottom row points in two rows of points that are set one above the 
other with respect to a choice of vertical. The strands are not individually 
knotted and they are disjoint from one another. See Figures 4, 5 and 6 for 
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illustrations of braids and moves on braids. Braids can be multiplied by at- 
taching the bottom row of one braid to the top row of the other braid. Taken 
up to ambient isotopy, fixing the endpoints, the braids form a group under this 
notion of multiplication. In Figure 4 we illustrate the form of the basic gener- 
ators of the braid group, and the form of the relations among these generators. 
Figure 5 illustrates how to close a braid by attaching the top strands to the 
bottom strands by a collection of parallel arcs. A key theorem of Alexander 
states that every knot or hnk can be represented as a closed braid. Thus the 
theory of braids is critical to the theory of knots and links. Figure 6 illustrates 
the famous Borromean Rings (a link of three unknotted loops such that any 
two of the loops are unlinked) as the closure of a braid. 

Let Bn denote the Artin braid group on n strands. We recall here that Bn 
is generated by elementary braids {si, • • • , with relations 

1. SiSj = SjSi for \i — j| > 1, 

2. SiSi+iSi = Si+iSiSi+i for i = 1, • • • n - 2. 

See Figure 4 for an illustration of the elementary braids and their relations. 
Note that the braid group has a diagrammatic topological interpretation, where 
a braid is an intertwining of strands that lead from one set of n points to 
another set of n points. The braid generators Sj are represented by diagrams 
where the i-th and (i + l)-th strands wind around one another by a single 
half-twist (the sense of this turn is shown in Figure 4) and all other strands 
drop straight to the bottom. Braids are diagrammed vertically as in Figure 4, 
and the products are taken in order from top to bottom. The product of two 
braid diagrams is accomplished by adjoining the top strands of one braid to 
the bottom strands of the other braid. 

In Figure 4 we have restricted the illustration to the four-stranded braid 
group B4. In that figure the three braid generators of B4 are shown, and then 
the inverse of the first generator is drawn. Following this, one sees the identities 
SiSi^ — 1 (where the identity element in B4 consists in four vertical strands), 
S1S2S1 — S2S1S2, and finally S1S3 = s^Si. 

Braids are a key structure in mathematics. It is not just that they are a 
collection of groups with a vivid topological interpretation. From the algebraic 
point of view the braid groups Bn are important extensions of the symmetric 
groups Sn- Recall that the symmetric group Sn of all permutations of n distinct 
objects has presentation as shown below. 

1. — 1 for i — 1, ■ ■ ■ n — 1, 
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2. SiSj = SjSi for |i — j| > 1, 

3. SiSi+iSi = Si+iSiSi+i for i = 1, • • • n - 2. 

Thus Sn is obtained from Bn by setting the square of each braiding generator 
equal to one. We have an exact sequence of groups 

1 > B^ > Sji ^ 1 

exhibiting the Artin Braid group as an extension of the symmetric group. 

In the next sections we shall show how representations of the Artin Braid 
group are rich enough to provide a dense set of transformations in the unitary 
groups. Thus the braid groups are, in principle fundamental to quantum 
computation and quantum information theory. 



2 Quantum Mechanics and Quantum Compu- 
tation 

We shall quickly indicate the basic principles of quantum mechanics. The 
quantum information context encapsulates a concise model of quantum theory: 

The initial state of a quantum process is a vector \v) in a complex vector 
space H. Measurement returns basis elements (3 of H with probability 

\{P\v)\V{v\v) 

where {v \w) — v^w with the conjugate transpose of v. A physical process oc- 
curs in steps \v) — > U \v) — \Uv) where U is a unitary linear transformation. 

Note that since {Uv\Uw) = {v \ WU\w) = {v \ w) — when U is unitary, it 
follows that probability is preserved in the course of a quantum process. 

One of the details for any specific quantum problem is the nature of the 
unitary evolution. This is specified by knowing appropriate information about 
the classical physics that supports the phenomena. This information is used 
to choose an appropriate Hamiltonian through which the unitary operator is 
constructed via a correspondence principle that replaces classical variables with 
appropriate quantum operators. (In the path integral approach one needs a 
Langrangian to construct the action on which the path integral is based.) One 
needs to know certain aspects of classical physics to solve any specific quantum 
problem. 
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A key concept in the quantum information viewpoint is the notion of the 
superposition of states. If a quantum system has two distinct states \v) and 
\w), then it has infinitely many states of the form a\v) + b\w) where a and b 
are complex numbers taken up to a common multiple. States are "really" in 
the projective space associated with H. There is only one superposition of a 
single state 1^) with itself. 

Dirac [21] introduced the "bra -(c)-ket" notation {A\B) — A^B for the 
inner product of complex vectors A, B & H. He also separated the parts of 
the bracket into the bra < A \ and the ket \B). Thus 

{A\B) = {A\ \B) 

In this interpretation, the ket \B) is identified with the vector B & H, while the 
bra < A I is regarded as the element dual to A in the dual space H*. The dual 
element to A corresponds to the conjugate transpose of the vector A, and 
the inner product is expressed in conventional language by the matrix product 
A'^B (which is a scalar since B is a column vector). Having separated the bra 
and the ket, Dirac can write the "ket-bra" = AB'^ . In conventional 

notation, the ket-bra is a matrix, not a scalar, and we have the following 
formula for the square of P = | : 

p2 = \A){B\\A){B I = A{B^A)B^ = {B^A)AB^ = {B \A)P. 

The standard example is a ket-bra P = \A){A\ where (A |A) = 1 so that 
p2 _ p fjj^gj^ P is a projection matrix, projecting to the subspace of H that 
is spanned by the vector \A). In fact, for any vector |P) we have 

P\B)^\A){A\\B)^\A){A\B)^{A\B)\A). 
If {|Ci), IC2), • • • |C„)} is an orthonormal basis for H, and 

Pi = I ) (Cj I , 

then for any vector \A) we have 

\A) = {CM)\Cl) + --- + {Cn\A)\Cn). 

Hence 

(P \A) = (P \A) + --- + {B \Cn){Cn \A) 

One wants the probability of starting in state \A) and ending in state |P). 
The probability for this event is equal to |(P This can be refined if we 

have more knowledge. If the intermediate states |Cj) are a complete set of 
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orthonormal alternatives then we can assume that {Ci \ Ci) = 1 for each i and 
that T,i\Ci){Ci\ = 1. This identity now corresponds to the fact that 1 is the 
sum of the probabihties of an arbitrary state being projected into one of these 
intermediate states. 

If there are intermediate states between the intermediate states this for- 
mulation can be continued until one is summing over all possible paths from 
A to B. This becomes the path integral expression for the amplitude {B\A). 

2.1 What is a Quantum Computer? 

A quantum, computer is, abstractly, a composition U of unitary transforma- 
tions, together with an initial state and a choice of measurement basis. One 
runs the computer by repeatedly initializing it, and then measuring the result 
of applying the unitary transformation U to the initial state. The results of 
these measurements arc then analyzed for the desired information that the 
computer was set to determine. The key to using the computer is the design 
of the initial state and the design of the composition of unitary transforma- 
tions. The reader should consult [65] for more specific examples of quantum 
algorithms. 

Let if be a given finite dimensional vector space over the complex numbers 
C. Let {Wo, Wi, Wn} be an orthonormal basis for H so that with \i) := \Wi) 
denoting Wi and {i\ denoting the conjugate transpose of we have 

where Sij denotes the Kronecker delta (equal to one when its indices are equal 
to one another, and equal to zero otherwise). Given a vector v in H let 
I^P := (^1^). Note that {i\v is the i-th coordinate of v. 

An measurement of v returns one of the coordinates \i) of v with probability 
K^lvp. This model of measurement is a simple instance of the situation with a 
quantum mechanical system that is in a mixed state until it is measured. The 
result of measurement is to put the system into one of the basis states. 

When the dimension of the space H is two (n = 1), a vector in the space 
is called a qubit. A qubit represents one quantum of binary information. On 
meausurement, one obtains either the ket |0) or the ket This constitutes 
the binary distinction that is inherent in a qubit. Note however that the 
information obtained is probabilistic. If the qubit is 

|V') = «|0)+/3 
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then the ket |0) is observed with probabihty lap, and the ket |1) is observed 
with probabihty In speaking of an ideahzed quantum computer, we do not 
specify the nature of measurement process beyond these probabihty postulates. 

In the case of general dimension n of the space H, we will call the vectors 
in H qudits. It is quite common to use spaces H that are tensor products 
of two-dimensional spaces (so that all computations are expressed in terms of 
qubits) but this is not necessary in principle. One can start with a given space, 
and later work out factorizations into qubit transformations. 

A quantum computation consists in the application of a unitary transfor- 
mation U to an initial qunit ip = ao|0) -|- ... -|- an\n) with \ip\'^ = 1, plus a 
measurement of Uip. A measurement of Uip returns the ket \i) with probabil- 
ity |(i|C/'0p. In particular, if we start the computer in the state then the 
probability that it will return the state \j) is |(j|C/|i)p. 

It is the necessity for writing a given computation in terms of unitary 
transformations, and the probabilistic nature of the result that characterizes 
quantum computation. Such computation could be carried out by an idealized 
quantum mechanical system. It is hoped that such systems can be physically 
realized. 

2.2 Universal Gates 

A two- qubit gate G is a unitary linear mapping G : V <^V — > V where V is 
a two complex dimensional vector space. We say that the gate G is universal 
for quantum computation (or just universal) if G together with local unitary 
transformations (unitary transformations from V to V) generates all unitary 
transformations of the complex vector space of dimension 2" to itself. It is well- 
known [65] that CNOT is a universal gate. (On the standard basis, GNOT is 
the identity when the first qubit is 0, and it flips the second qbit, leaving the 
first alone, when the first qubit is 1.) 

A gate G, as above, is said to be entangling if there is a vector 

\al3) = \a)®\l3) ^V®V 

such that G\a(3) is not decomposable as a tensor product of two qubits. Under 
these circumstances, one says that G\al3) is entangled. 

In [16], the Brylinskis give a general criterion of G to be universal. They prove 
that a two- qubit gate G is universal if and only if it is entangling. 
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Remark. A two-qubit pure state 



10) = a|00) + 6|01) + c|10) + d\n) 



is entangled exactly when {ad — be) ^ 0. It is easy to use this fact to check 
when a specific matrix is, or is not, entangling. 

Remcirk. There are many gates other than CNOT that can be used as 
universal gates in the presence of local unitary transformations. Some of these 
are themselves topological (unitary solutions to the Yang-Baxter equation, see 
[52, 81]) and themselves generate representations of the Artin Braid Group. 
Replacing CNOT by a solution to the Yang-Baxter equation does not place 
the local unitary transformations as part of the corresponding representation 
of the braid group. Thus such substitutions connote only a partial solution 
to creating topological quantum computation. In this paper we are concerned 
with braid group representations that include all aspects of the unitary group. 
Accordingly, in the next section we shall first examine, how the braid group 
on three strands can be represented as local unitary transformations. 



3 SU{2) Representations of the Artin Braid 
Group 



The purpose of this section is to determine all the representations of the three 
strand Artin braid group to the special unitary group SU (2) and concomi- 
tantly to the unitary group U{2). One regards the groups SU{2) and U{2) as 
acting on a single qubit, and so U{2) is usually regarded as the group of local 
unitary transformations in a quantum information setting. If one is looking 
for a coherent way to represent all unitary transformations by way of braids, 
then U{2) is the place to start. Here we will show that there are many rep- 
resentations of the three-strand braid group that generate a dense subset of 
U{2). Thus it is a fact that local unitary transformations can be "generated 
by braids" in many ways. 

We begin with the structure of SU{2). A matrix in SU{2) has the form 



where z and w are complex numbers, and z denotes the complex conjugate of 
z. To be in SU{2) it is required that Det(M) = 1 and that = M'^ where 
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Det denotes determinant, and M'^ is the conjugate transpose of M. Thus if 
z — a + bi and w = c + di where a, 6, c, d are real numbers, and = — 1, then 



_ f a + bi c + di \ 
\ —c + di a — hi j 

with + 6^ + + = 1. It is convenient to write 

and to abbreviate this decomposition as 

M = a + hi + cj + dk 



where 



1 i I -i i '-^ I -1 i I i 



so that 
and 



i = j = k = ijk = —1 

ij = kjk = i,ki = j 

ji = -k, kj = -i, ik = -j. 

The algebra of l,i,j,k is called the quaternions after William Rowan Hamil- 
ton who discovered this algebra prior to the discovery of matrix algebra. Thus 
the unit quaternions are identified with SU{2) in this way. We shall use this 
identification, and some facts about the quaternions to find the SU (2) repre- 
sentations of braiding. First we recall some facts about the quaternions. 



1. Note that ii q = a + hi + cj + dk (as above), then q^ — a — bi — cj — dk 
so that gg'^ ^ a"^ + b'^ + + d^ ^ 1. 

2. A general quaternion has the form q — a + bi + cj + dk where the value 
of qq'^ — + b'^ + + d'^ , is not fixed to unity. The length of q is by 
definition ^Jqq^. 

3. A quaternion of the form ri + sj + tk for real numbers r, s, t is said to 
be a pure quaternion. We identify the set of pure quaternions with the 
vector space of triples (r, s, t) of real numbers R^. 
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4. Thus a general quaternion has the form q = a -\- bu where -u is a pure 
quaternion of unit length and a and b are arbitrary real numbers. A unit 
quaternion (element of SU{2)) has the addition property that = 1. 

5. If u is a pure unit length quaternion, then -u^ = —1. Note that the 
set of pure unit quaternions forms the two-dimensional sphere 5*^ — 
{(r,s,t)|r2 + s2 + i2 = l} in 

6. U u,v are pure quaternions, then 

uv — —u • V -\-ux V 

whre u ■ V is the dot product of the vectors u and and u x v is the 
vector cross product of u and v. In fact, one can take the definition of 
quaternion multiplication as 

(a + bu){c + dv) = ac + bc{u) + ad{y) + bd{—u ■ v + u x v), 

and all the above properties are consequences of this definition. Note 
that quaternion multiplication is associative. 

7. Let g — a + bu he a, unit length quaternion so that v? — —1 and a — 
cos{9/2),b = sin{9/2) for a chosen angle 9. Define (f)g : — > R^ by 
the equation (l>g{P) — gPg\ for P any point in i?^, regarded as a pure 
quaternion. Then (j)g is an orientation preserving rotation of R^ (hence 
an element of the rotation group 5*0(3)). Specifically, 0g is a rotation 
about the axis u by the angle 9. The mapping 

(j) : SU{2) — > S0{2,) 

is a two-to-one surjective map from the special unitary group to the 
rotation group. In quaternionic form, this result was proved by Hamilton 
and by Rodrigues in the middle of the nineteeth century. The specific 
formula for 4>g{P) as shown below: 

(f)g{P) = gPg-^ = (a^ - b^)P + 2ab{P x w) + 2(P • uyu. 

We want a representation of the three-strand braid group in SU{2). This 
means that we want a homomorphism p : — > SU{2), and hence we want 
elements g = p{si) and h = p{s2) in SU{2) representing the braid group 
generators Si and s^- Since S1S2S1 = S2S1S2 is the generating relation for i?3, 
the only requirement on g and h is that ghg — hgh. We rewrite this relation 
as h~^gh — ghg~^, and analyze its meaning in the unit quaternions. 
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Suppose that g = a + bu and h = c + dv where u and v are unit pure 
quaternions so that + 6^ = 1 and + dP = 1. then ghg~^ = c + d(f)g{v) and 
h'^gh — a + h(f)h-i{u). Thus it follows from the braiding relation that a — 
b = ±d, and that (pg{v) — ±(f)fi-i{u). However, in the case where there is a 
minus sign we have g = a + bu and h — a — bv — a-\- b{—v). Thus we can now 
prove the following Theorem. 

Theorem 1. If g — a + bu and h — c + dv are pure unit quaternions, then, 
without loss of generality, the braid relation ghg — hgh is true if and only if 

h — a + bv, and (f)g{v) = (j)h-i{u). Furthermore, given that g = a + bu and 
h = a + bv, the condition (j)g{v) = 0^-1 (m) is satisfied if and only ifu-v = "^~2 
when u ^ V. If u = V then then g = h and the braid relation is trivially 
satisfied. 

Proof. We have proved the first sentence of the Theorem in the discussion 
prior to its statement. Therefore assume that g = a + bu, h = a + bv , and 
(f>g{v) = (ph-i (m). We have already stated the formula for (f)g{v) in the discussion 
about quaternions: 

(j)g{v) = gvg~^ = (a^ - b^)v + 2ab{v x u) + 2{v ■ u)b'^u. 

By the same token, we have 

(t)h-i{u) = h'^uh = (a^ - b^)u + 2ab{u x -v) + 2{u ■ {-v))b'^{-v) 

= (a^ - b'^)u + 2ab{v x u) + 2{v ■ u)b'^{v). 
Hence we require that 

(a^ - b'^)v + 2{v ■ u)b'^u = (a^ - b'^)u + 2{v ■ u)b'^{v). 

This equation is equivalent to 

2{u ■ v)b'^{u -v) = (a^ - b'^){u - v). 

li u ^ V, then this implies that 

a^-b^ 

u ■ V — — — . 

262 

This completes the proof of the Theorem. □ 
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An Example. Let 



9 = 



a + bi 



where a = cos{9) and h = sin{6). Let 



s^)i + 2csk] 



where c^-\-s^ — 1 and (? — s^ — ^ . Then we can reexpress g and h in matrix 
form as the matrices G and H. Instead of writing the exphcit form of if, we 
write H — FGF^ where F is an element of SU (2) as shown below. 



This representation of braiding where one generator G is a simple matrix of 
phases, while the other generator H = FGF^ is derived from G by conjugation 

by a unitary matrix, has the possibility for generalization to representations of 
braid groups (on greater than three strands) to SU (n) or U (n) for n greater 
than 2. In fact we shall see just such representations constructed later in this 
paper, by using a version of topological quantum field theory. The simplest 
example is given by 



where r^+r = 1. Then g and h satisfy ghg — hgh and generate a representation 
of the three-strand braid group that is dense in SU{2). We shall call this the 
Fibonacci representation of B3 to SU{2). 

Density. Consider representations of into SU{2) produced by the method 
of this section. That is consider the subgroup SU[G, H] of SU{2) generated by 
a pair of elements {(?, h} such that ghg = hgh. Wc wish to understand when 
such a representation will be dense in SU{2). We need the following lemma. 

Lemma. e"*e''-'e'^* = cos{h)e''^"'^^^ + sin{h)e^^"'~'^'> j. Hence any element ofSU{2) 
can be written in the form e'^'e^'^e'^''. for appropriate choices of angles a, b, c. In 
fact, if u and v are linearly independent unit vectors in R^, then any element 
of SU (2) can be written in the form 



for appropriate choices of the real numbers a, b, c. 
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Proof. It is easy to check that 

This completes the verification of the identity in the statement of the Lemma. 
Let V be any unit direction in and A an arbitrary angle. We have 

e"'^ = cos (A) + sin{X)v, 

and 

V — r + si + {p + qi)j 
where + s^ + + — 1. So 

e^^ = cos(A) + sin{X)[r + si] + sin{X)\p + qi]j 

— [(cos(A) + sm(A)r) + sin{X)si] + [sin{X)p + sin{X)qi]j. 

By the identity just proved, we can choose angles a, b, c so that 

Hence 

cos (6) 6^^"+''^ = (cos(A) + sm(A)r) + sin{X)si 

and 

sin(b)e'^"'~''^ = sin{X)p + sin{X)qi. 

Suppose we keep v fixed and vary A. Then the last equations show that this 
will result in a full variation of b. 



Now consider 



gia'giJAgic' _ ^ia' ^ia^jb^ic^ib' _ ^i{a'+a) ^jb^i{c+c') 



By the basic identity, this shows that any element of SU (2) can be written in 
the form 



Then, by applying a rotation, we finally conclude that if u and v are linearly 
independent unit vectors in R^, then any element of SU (2) can be written in 
the form 



for appropriate choices of the real numbers a,b,c. □ 
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This Lemma can be used to verify density of a representation, by finding 
two elements A and B in tlie representation sucli tliat tlie powers of A are dense 
in the rotations about its axis, and the powers of B are dense in the rotations 
about its axis, and such that the axes of A and B are hnearly independent in 
R^. Then by the Lemma the set of elements A°'^'^B^A°'^'^ are dense in SU{2). 
It follows for example, that the Fibonacci representation described above is 
dense in SU{2), and indeed the generic representation of B^ into SU{2) will 
be dense in SU{2). Our next task is to describe representations of the higher 
braid groups that will extend some of these unitary repressentations of the 
three-strand braid group. For this we need more topology. 



4 The Bracket Polynomial and the Jones Poly- 
nomial 

We now discuss the Jones polynomial. We shall construct the Jones polynomial 
by using the bracket state summation model [34]. The bracket polynomial, 
invariant under Reidmeister moves II and III, can be normalized to give an 

invariant of all three Reidcmcistcr moves. This normalized invariant, with a 
change of variable, is the Jones polynomial [32, 33]. The Jones polynomial was 
originally discovered by a different method than the one given here. 

The bracket polynomial , < K > = < K > (A) , assigns to each unoriented 
link diagram K a Laurent polynomial in the variable A, such that 

1. If and K' are regularly isotopic diagrams, then < K > = < K' >. 

2. U KUO denotes the disjoint union of K with an extra unknotted and un- 
linked component O (also called 'loop' or 'simple closed curve' or 'Jordan 
curve'), then 

< KUO >^5 < K >, 

where 

5^-A^- A-\ 

3. < K > satisfies the following formulas 

<X>=^<>^> +^"^ <)(> 
<X> = ^"'<->+^ <)(>, 



21 



where the small diagrams represent parts of larger diagrams that are identical 
except at the site indicated in the bracket. We take the convention that the 
letter chi, denotes a crossing where the curved line is crossing over the 
straight segment. The barred letter denotes the switch of this crossing, where 
the curved line is undercrossing the straight segment. See Figure 7 for a graphic 
illustration of this relation, and an indication of the convention for choosing 
the labels A and A^^ at a given crossing. 





Figure 7 - Bracket Smoothings 



It is easy to see that Properties 2 and 3 define the calculation of the bracket 

on arbitrary link diagrams. The choices of coefficients (^4 and A~^) and the 
value of d make the bracket invariant under the Reidemeister moves II and III. 
Thus Property 1 is a consequence of the other two properties. 

In computing the bracket, one finds the following behaviour under Reide- 
meister move I: 

< 7 >= -A^ <^> 

and 

< 7 >= -A~^ <^> 

where 7 denotes a curl of positive type as indicated in Figure 8, and 7 indicates 
a curl of negative type, as also seen in this figure. The type of a curl is the sign 
of the crossing when we orient it locally. Our convention of signs is also given 
in Figure 8. Note that the type of a curl does not depend on the orientation 
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we choose. The small arcs on the right hand side of these formulas indicate 
the removal of the curl from the corresponding diagram. 

The bracket is invariant under regular isotopy and can be normalized to an 
invariant of ambient isotopy by the definition 

where we chose an orientation for K, and where w{K) is the sum of the crossing 
signs of the oriented link K. w{K) is called the writhe of K. The convention 
for crossing signs is shown in Figure 8. 




Figure 8 - Crossing Signs and Curls 

One useful consequence of these formulas is the following switching formula 

A<x> <X>= (^^ - ^~^) <^> ■ 

Note that in these conventions the 74-smoothing of % is while the A- 
smoothing of % is )(. Properly interpreted, the switching formula above says 
that you can switch a crossing and smooth it either way and obtain a three 
diagram relation. This is useful since some computations will simplify quite 
quickly with the proper choices of switching and smoothing. Remember that 
it is necessary to keep track of the diagrams up to regular isotopy (the equiv- 
alence relation generated by the second and third Reidemeister moves). Here 
is an example. View Figure 9. 
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Figure 9 — Trefoil and Two Relatives 

Figure 9 shows a trefoil diagram K, an unknot diagram U and another unknot 
diagram U' . Applying the switching formula, we have 

A-^ <K > -A<U >^ - A^) <U'> 

and < [/ >= -A^ and < U' >= (-A'^y = A'^. Thus 

A-^ <K> -A{-A^) = {A-^ - A^)A-^ 

Hence 

A-^ <K>^ -A^ + A-^ - A-^. 

Thus 

<K>= -A^ - A-^ + A-\ 
This is the bracket polynomial of the trefoil diagram K. 

Since the trefoil diagram K has writhe w{K) — 3, we have the normalized 
polynomial 

/^(A) = <K>= -A-\-J^' - A-'^ + A-^) = A-" + A-^^ - A~^\ 

The bracket model for the Jones polynomial is quite useful both theoreti- 
cally and in terms of practical computations. One of the neatest applications 
is to simply compute, as wc have done, fK{A) for the trefoil knot K and de- 
termine that fxi^) is not equal to /^(A^^) = f_K{A). This shows that the 
trefoil is not ambient isotopic to its mirror image, a fact that is much harder 
to prove by classical methods. 
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The State Summation. In order to obtain a closed formula for tlie bracket, 
we now describe it as a state summation. Let K be any unoriented link 
diagram. Define a state, S, of X to be a choice of smoothing for each crossing 
of K. There are two choices for smoothing a given crossing, and thus there are 
2^ states of a diagram with crossings. In a state we label each smoothing 
with A or A^^ according to the left-right convention discussed in Property 3 
(see Figure 7). The label is called a vertex weight of the state. There are 
two evaluations related to a state. The first one is the product of the vertex 
weights, denoted 

< K\S > . 

The second evaluation is the number of loops in the state denoted 

ll-^ll- 

Define the state summation, < >, by the formula 

<K> = Y,<K\S>5\\^\\-\ 

s 

It follows from this definition that < K > satisfies the equations 

<X> = ^<>^ >+^"' <)(>, 
<KU0> = 5 <K >, 

< O > = 1. 

The first equation expresses the fact that the entire set of states of a given 
diagram is the union, with respect to a given crossing, of those states with 
an A-type smoothing and those with an A^^-type smoothing at that crossing. 
The second and the third equation are clear from the formula defining the state 
summation. Hence this state summation produces the bracket polynomial as 
we have described it at the beginning of the section. 

Remcirk. By a change of variables one obtains the original Jones polynomial, 
VK{t)i for oriented knots and links from the normalized bracket: 

VKit) = fK{t--^). 

Remark. The bracket polynomial provides a connection between knot theory 
and physics, in that the state summation expression for it exhibits it as a 
generalized partition function defined on the knot diagram. Partition functions 
are ubiquitous in statistical mechanics, where they express the summation 
over all states of the physical system of probability weighting functions for the 
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individual states. Such physical partition functions contain large amounts of 
information about the corresponding physical system. Some of this information 
is directly present in the properties of the function, such as the location of 
critical points and phase transition. Some of the information can be obtained 
by differentiating the partition function, or performing other mathematical 
operations on it. 

There is much more in this connection with statistical mechanics in that 

the local weights in a partition function are often expressed in terms of solu- 
tions to a matrix equation called the Yang-Baxter equation, that turns out to 
fit perfectly invariance under the third Reidemeister move. As a result, there 
are many ways to define partition functions of knot diagrams that give rise to 
invariants of knots and links. The subject is intertwined with the algebraic 
structure of Hopf algebras and quantum groups, useful for producing system- 
atic solutions to the Yang-Baxter equation. In fact Hopf algebras are deeply 
connected with the problem of constructing invariants of three-dimensional 
manifolds in relation to invariants of knots. We have chosen, in this survey 
paper, to not discuss the details of these approaches, but rather to proceed 
to Vassiliev invariants and the relationships with Wittcn's functional integral. 
The reader is referred to [34, 35, 36, 37, 40, 41, 3, 32, 33, 42, 69, 70, 74, 75] for 
more information about relationships of knot theory with statistical mechan- 
ics, Hopf algebras and quantum groups. For topology, the key point is that 
Lie algebras can be used to construct invariants of knots and links. 



4.1 Quantum Computation of the Jones Polynomial 

Quantum algorithms for computing the Jones polynomial have been discussed 
elsewhere. See [46, 52, 1, 54, 2, 79]. Here, as an example, we give a local 
unitary representation that can be used to compute the Jones polynomial for 
closures of 3-braids. We analyse this representation by making explicit how 
the bracket polynomial is computed from it, and showing how the quantum 
computation devolves to finding the trace of a unitary transformation. 

The idea behind the construction of this representation depends upon the 
algebra generated by two single qubit density matrices (ket-bras). Let \v) 
and \w) be two qubits in V, a complex vector space of dimension two over 
the complex numbers. Let P — \v){v\ and Q — \w){w\ be the corresponding 
ket-bras. Note that 

<5' = kr<5, 
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PQP = \ {v\w)\^P, 
QPQ ^ \{v\w)\^Q. 

P and Q generate a representation of the Temperlcy-Lieb algebra (See Section 
5 of the present paper). One can adjust parameters to make a representation 
of the three-strand braid group in the form 

si I — > rP + si, 

S2 I — > tQ + ul, 

where / is the identity mapping on V and r, s, t, u are suitably chosen scalars. 
In the following we use this method to adjust such a representation so that it 
is unitary. Note also that this is a local unitary representation of S3 to U{2). 
We leave it as an exersise for the reader to verify that it fits into our general 
classification of such representations as given in section 3 of the present paper. 

The representation depends on two symmetric but non- unitary matrices Ui 
and U2 with 



d 




= d\w){w\ 



and 



d- 



d- 



d-' 



d\v){v\ 



where w — (1, 0), and v = {d ^, \/l — d~^), assuming the entries of v are real. 
Note that Uf = dUi and f/| = dUi. Moreover, U1U2U1 = Ui and U2U1U2 = Ui. 
This is an example of a specific representation of the Temperlcy-Lieb algebra 
[34, 46]. The desired representation of the Artin braid group is given on the 
two braid generators for the three strand braid group by the equations: 

^82) ^ AI + A-^U2. 
Here / denotes the 2x2 identity matrix. 

For any A with d = —A"^ — A~'^ these formulas define a representation of the 
braid group. With A = e*^, we have d = —2cos{29) . We find a specific range 
of angles |^| < 7r/6 and — 7r| < 7r/6 that give unitary representations of the 
three-strand braid group. Thus a specialization of a more general represention 
of the braid group gives rise to a continuous family of unitary representations 
of the braid group. 
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Note that the traces of these matrices are given by the formulas tr{Ui) = 
tr{U2) = d while tr{UiU2) = tr{U2Ui) = 1. If 6 is any braid, let I{h) denote the 
sum of the exponents in the braid word that expresses h. For h a three-strand 
braid, it follows that 

$(6) = + n(6) 

where / is the 2x2 identity matrix and 11 (fe) is a sum of products in the 

Temperley-Lieb algebra involving Ui and U2- Since the Temperley-Lieb algebra 
in this dimension is generated by /,f/i, U2-, U1U2 and U2U1, it follows that the 
value of the bracket polynomial of the closure of the braid b, denoted < 6 > , 
can be calculated directly from the trace of this representation, except for the 
part involving the identity matrix. The bracket polynomial evaluation depends 
upon the loop counts in the states of the closure of the braid, and these loop 
counts correspond to the traces of the non-identity Temperley-Lieb elements. 
Note that the closure of the three-strand diagram for the identity braid in B3 
has bracket polynomial d^. The result is the equation 

<b>^A^^''^d'' + tr{U{b)) 

where b denotes the standard braid closure of 6, and the sharp brackets denote 
the bracket polynomial. Since the trace of the 2x2 identity matrix is 2, we 
see that 

<b>^trmb))+A'^^\d^ -2). 

It follows from this calculation that the question of computing the bracket 
polynomial for the closure of the three-strand braid b is mathematically equiv- 
alent to the problem of computing the trace of the unitary matrix $(6). 

The Hadamard Test 

In order to (quantum) compute the trace of a unitary matrix U, one can use 
the Hadamard test to obtain the diagonal matrix elements {ip\U\ip) of U. The 
trace is then the sum of these matrix elements as 1-0) runs over an orthonormal 
basis for the vector space. We first obtain 

as an expectation by applying the Hadamard gate H 

H\o) = ^m + \^)) 
H\i) = - 
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to the first qubit of 



Cuo{H(^ 1)\0M = -^(|0) ® 1^) + |1) ® ui^p). 

Here Cu denotes controlled U, acting as U when the control bit is |1) and the 
identity mapping when the control bit is |0). We measure the expectation for 
the first qubit |0) of the resulting state 

^{H\o) iv^) +H\i)^ um = ^((|o) + |i)) ® \^) + (|o) - |i)) ® um 
= Im ® m + um + \i) ® (iv'> - u\m- 

This expectation is 

The imaginary part is obtained by applying the same procedure to 

i=(|o)0|v^)-z|i)0t/|V^) 

This is the method used in [1], and the reader may wish to contemplate its 
efficiency in the context of this simple model. Note that the Hadamard test 
enables this quantum computation to estimate the trace of any unitary ma- 
trix U by repeated trials that estimate individual matrix entries {t/j\U\t/j). We 
shall return to quantum algorithms for the Jones polynomial and other knot 
polynomials in a subsequent paper. 



5 Quantum Topology, Cobordism Categories, 
Temperley-Lieb Algebra and Topological Quan- 
tum Field Theory 

The purpose of this section is to discuss the general idea behind topological 
quantum field theory, and to illustrate its application to basic quantum me- 
chanics and quantum mechanical formalism. It is useful in this regard to have 
available the concept of category, and we shall begin the section by discussing 
this far-reaching mathematical concept. 
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Definition. A category Cat consists in two related collections: 

1. Obj{Cat), the objects of Cat, and 

2. Morph{Cat), the morphisms of Cat. 
satisfying the following axioms: 

1. Each morphism / is associated to two objects of Cat, the domain of f 
and the codomain of f. Letting A denote the domain of / and B denote 
the codomain of /, it is customary to denote the morphism / by the 
arrow notation / : A — > B. 

2. Given / : A — > B and g : B — > C where A, B and C are objects of 
Cat, then there exists an associated morphism go f : A — > C called the 
composition of / and g. 

3. To each object A of Cat there is a unique identity morphism Ij^ : A — > A 
such that 1a o/ = / for any morphism / with codomain A, and goljt^ = g 
for any morphism g with domain A. 

4. Given three morphisms / : A — > B, g : B — > C and h : C — > D, then 
composition is associative. That is 

{ho g) o f = ho {g o f). 

If Cati and Cat2 are two categories, then a functor F : Cati — > Cat2 consists 
in functions Fq : Obj{Cati) — > Obj{Cat2) and Fm : Morph{Cati) — > 
Morph{Cat2) such that identity morphisms and composition of morphisms 
are preserved under these mappings. That is (writing just F for Fq and Fm), 

1. F{Ia) = 1f(A), 

2. F{f:A^B)^ F{f) : F{A) F{B), 

3. F{gof) = F{g)oF{f). 

A functor F : Cati — > Cat2 is a structure preserving mapping from one 
category to another. It is often convenient to think of the image of the functor 
F as an interpretation of the first category in terms of the second. We shall 
use this terminology below and sometimes refer to an interpretation without 
specifying all the details of the functor that describes it. 
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The notion of category is a broad mathematical concept, encompassing 
many fields of mathematics. Thus one has the category of sets where the 
objects are sets (collections) and the morphisms are mappings between sets. 
One has the category of topological spaces where the objects are spaces and 
the morphisms are continuous mappings of topological spaces. One has the 
category of groups where the objects arc groups and the morphisms arc homo- 
morphisms of groups. Functors are structure preserving mappings from one 
category to another. For example, the fundamental group is a functor from 
the category of topological spaces with base point, to the category of groups. 
In all the examples mentioned so far, the morphisms in the category are re- 
strictions of mappings in the category of sets, but this is not necessarily the 
case. For example, any group G can be regarded as a category, Cat{G), with 
one object *. The morphisms from * to itself are the elements of the group 
and composition is group multiplication. In this example, the object has no 
internal structure and all the complexity of the category is in the morphisms. 

The Artin braid group -B„ can be regarded as a category whose single 
object is an ordered row of points [n] — {1, 2, 3, n}. The morphisms are the 
braids themselves and composition is the multiplication of the braids. The 
ordered row of points is interpreted as the starting and ending row of points 
at the bottom and the top of the braid. In the case of the braid category, 
the morphisms have both external and internal structure. Each morphism 
produces a permutation of the ordered row of points (corresponding to the 
begiinning and ending points of the individual braid strands) , and weaving of 
the braid is extra structure beyond the object that is its domain and codomain. 
Finally, for this example, we can take all the braid groups {n a positive 
integer) under the wing of a single category, Cat{B), whose objects are all 
ordered rows of points [n] , and whose morphisms are of the form b : [n] — > [n] 
where & is a braid in B^. The reader may wish to have morphisms between 
objects with different n. We will have this shortly in the Temper ley- Lieb 
category and in the category of tangles. 

The n-Cobordism Category, Cob[n\, has as its objects smooth manifolds of 
dimension n, and as its morphisms, smooth manifolds M"+^ of dimension n + 1 
with a partition of the boundary, (9M"+^, into two collections of n-manifolds 
that we denote by L{M"'^^) and i?(M""'"^). We regard M^''^^ as a morphism 
from L(M"+i) to R{M''+^) 

M"+^ : L(M'*+^) — > i?(M"+i). 

As we shall see, these cobordism categories are highly significant for quantum 
mechanics, and the simplest one, Co6[0] is directly related to the Dirac notation 
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of bras and kets and to the Temperley-Lieb algebara. We shall concentrate 
in this section on these cobordism categories, and their relationships with 
quantum mechanics. 

One can choose to consider either oriented or non-oriented manifolds, and 
within unoriented manifolds there are those that are orientable and those that 
are not orientable. In this section we will implicitly discuss only orientable 
manifolds, but we shall not specify an orientation. In the next section, with 
the standard definition of topological quantum field theory, the manifolds will 
be oriented. The definitions of the cobordism categories for oriented manifolds 
go over mutatis mutandis. 

Lets begin with Co6[0]. Zero dimensional manifolds are just collections 
of points. The simplest zero dimensional manifold is a single point p. We 
take p to be an object of this category and also *, where * denotes the empty 
manifold (i.e. the empty set in the category of manifolds). The object * occurs 
in Coh[n] for every n, since it is possible that either the left set or the right set 
of a morphism is empty. A line segment S with boundary points p and 5 is a 
morphism from pto q. 

S:p^q 

See Figure 10. In this figure we have illustrated the morphism from p to p. 
The simplest convention for this category is to take this morphism to be the 
identity. Thus if we look at the subcategory of Co6[0] whose only object is p, 
then the only morphism is the identity morphism. Two points occur as the 
boundary of an interval. The reader will note that Co6[0] and the usual arrow 
notation for morphisms are very closely related. This is a place where notation 
and mathematical structure share common elements. In general the objects of 
Co6[0] consist in the empty object * and non-empty rows of points, symbolized 

by 

p(g)p(g)---(g)p(g)p. 

Figure 10 also contains a morphism 

p <8)p — > * 

and the morphism 

* — > p®p. 

The first represents a cobordism of two points to the empty set (via the bound- 
ing curved interval). The second represents a cobordism from the empty set 
to two points. 
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• • -^-^ Identity 

; :>-^ p(g)p 

czi:^^ * >P(8)p 

Figure 10 - Elementary Cobordisms 

In Figure 11, we have indicated more morphisms in Co6[0], and we have 
named the morphisms just discussed as 

|0) :p(g)p — > *, 

(©I : * — > p®p. 

The point to notice is that the usual conventions for handling Dirac bra-kets 
arc essentially the same as the compostion rules in this topological category. 
Thus in Figure 11 we have that 

(e| o \n) = (e|fi) : * — >* 

represents a cobordism from the empty manifold to itself. This cobordism is 
topologically a circle and, in the Dirac formalism is interpreted as a scalar. 
In order to interpret the notion of scalar we would have to map the cobor- 
dism category to the category of vector spaces and linear mappings. Wc shall 
discuss this after describing the similarities with quantum mechanical formal- 
ism. Nevertheless, the reader should note that if V is a vector space over the 
complex numbers C, then a linear mapping from C to C is determined by the 
image of 1, and hence is characterized by the scalar that is the image of 1. In 
this sense a mapping C — > C can be regarded as a possible image in vector 
spaces of the abstract structure (0|f^) : * — > *■ It is therefore assumed that 
in Co6[0] the composition with the morphism commutes with any other 

morphism. In that way (0|f^) behaves like a scalar in the cobordism category. 
In general, an 1 manifold without boundary behaves as a scalar in Cob[n], 
and if a manifold M""*"^ can be written as a union of two submanifolds L""*"^ 
and so that that an n-manifold is their common boundary: 
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with 



then, we can write 



and (M"+^) will be a scalar (morphism that commutes with all other mor- 
phisms) in the category Co6[n]. 



» • Identity 

czi: ^ <0i 

c ^ ^ <0|^1> 

= U 



• — 




=5c=: 


• — 
• — 






• — 




=>c=:: 


c 







u u = i^^> <e\Q.> <0i 

= <0|^l> 1^1x01 

= <0|^i> u 



Figure 11 - Bras, Kets and Projectors 
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su = us = u 

Figure 12 - Permutations 




{l^^> <0i}(g) 1 = p 
10{l^^><0l} = Q 



^ — :> -^{i(g)<0|} {l£l>0 1} 



Figure 13 - Projectors in Tensor Lines and Elementciry Topology 

Getting back to the contents of Figure 11, note how the zero dimensional 
cobordism category has structural parallels to the Dirac ket-bra formalism 

u^\n){e\ 

uu = |Q)(e|Q)(e| = (e|fi)|Q)(e| = {e\n)u. 

In the cobordism category, the bra-ket and ket-bra formalism is seen as pat- 
terns of connection of the one-manifolds that realize the cobordisms. 
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Now view Figure 12. This Figure illustrates a morphism S in Co6[0] that 
requires two crossed line segments for its planar representation. Thus S can 
be regarded as a non-trivial permutation, and — I where / denotes the 
identity morphisms for a two-point row. Prom this example, it is clear that 
Co6[0] contains the structure of all the symmetric groups and more. In fact, if 
we take the subcateogry of Co6[0] consisting of all morphisms from [n] to [n] 
for a fixed positive integer n, then this gives the well-known Brauer algebra 
(see [12]) extending the symmetric group by allowing any connections among 
the points in the two rows. In this sense, one could call Co&[0] the Brauer 
category. We shall return to this point of view later. 

In this section, we shall be concentrating on the part of Co6[0] that does not 
involve permutations. This part can be characterized by those morphisms that 
can be represented by planar diagrams without crosssings between any of the 
line segments (the one- manifolds). We shall call this crossingless subcategory 
of Co6[0] the Temperley-Lieb Category and denote it by CatTL. In CatTL we 
have the subcategory TL[n] whose only objects are the row of n points and the 
empty object *, and whose morphisms can all be represented by configurations 
that embed in the plane as in the morphisms P and Q in Figure 13. Note that 
with the empty object *, the morphism whose diagram is a single loop appears 
in TL[n] and is taken to commute with all other morphisms. 

The Temperley-Lieb Algebra, AlgTL[n\ is generated by the morphisms in 
TL[n] that go from [n] to itself. Up to multiplication by the loop, the product 
(composition) of two such morphisms is another flat morphism from [n] to 
itself. For algebraic purposes the loop * — > * is taken to be a scalar algebraic 
variable 6 that commutes with all elements in the algebra. Thus the equation 

uu = {e\n)u. 

becomes 

UU = 5U 

in the algebra. In the algebra we are allowed to add morphisms formally and 
this addition is taken to be commutative. Initially the algebra is taken with 
coefficients in the integers, but a different commutative ring of coefficients can 
be chosen and the value of the loop may be taken in this ring. For example, 
for quantum mechanical applications it is natural to work over the complex 
numbers. The multiplicative structure of AlgTL[n] can be described by gen- 
erators and relations as follows: Let denote the identity morphism from [n] 
to [n]. Let Ui denote the morphism from [n] to [n] that connects k with k for 
k < i and k > i + l from one row to the other, and connects i to i -|- 1 in each 
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row. Then the algebra AlgTL[n] is generated by {/„, Ui, U2, ■ ■ ■ , Un-i} with 
relations 

UiUi+iUi = Ui 

UiUj = UjU, : \i-j\>l. 

These relations are illustrated for three strands in Figure 13. We leave the 
commuting relation for the reader to draw in the case where n is four or 
greater. For a proof that these are indeed all the relations, see [49] . 

Figures 13 and 14 indicate how the zero dimensional cobordism category 
contains structure that goes well beyond the usual Dirac formalism. By ten- 
soring the kct-bra on one side or another by identity morphisms, wc obtain 
the beginnings of the Temperley-Lieb algebra and the Tcmperley-Licb cate- 
gory. Thus Figure 14 ilhistratcs the morphisms P and Q obtained by such 
tensoring, and the relation PQP = P which is the same as U1U2U1 = Ui 

Note the composition at the bottom of the Figure 14. Here we see a com- 
position of the identity tensored with a ket, followed by a bra tensored with 
the identity. The diagrammatic for this association involves "straightening" 
the curved structure of the morphism to a straight line. In Figure 15 we have 
elaborated this situation even further, pointing out that in this category each 
of the morphisms (0| and 1^2) can be seen, by straightening, as mappings 
from the generating object to itself. We have denoted these corresponding 
morphisms by 6 and fl respectively. In this way there is a correspondence 
between morphisms p (g) p — > * and morphims p — > p. 

In Figure 15 we have illustrated the generalization of the straightening 
procedure of Figure 14. In Figure 14 the straightening occurs because the 
connection structure in the morphism of Co6[0] does not depend on the wan- 
dering of curves in diagrams for the morphisms in that category. Nevertheless, 
one can envisage a more complex interpretation of the morphisms where each 
one-manifold (line segment) has a label, and a multiplicity of morphisms can 
correspond to a single line segment. This is exactly what we expect in inter- 
pretations. For example, we can interpret the line segment [1] — > [1] as a 
mapping from a vector space V to itself. Then [1] — > [1] is the diagrammatic 
abstraction for V — > V, and there are many instances of linear mappings from 
V to V. 

At the vector space level there is a duality between mappings V<S>V — > C 
and linear maps V — > V. Specifically, let 

{|o),---,|m)} 
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be a basis for V. Then : V — > V is determined by 

(where we have used the Einstein summation convention on the repeated index 
j) corresponds to the bra 

{e\:V0V — >C 

defined by 

(0|zj) = 0,,. 

Given (0| : V — > C, we associate : V — > V in this way. 

Comparing with the diagrammatic for the category Co6[0], we say that 
: V — > V is obtained by straightening the mapping 

{&\:V®V — >C. 

Note that in this interpretation, the bras and kets are defined relative to the 
tensor product of V with itself and [2] is interpreted as V (S)V. If we interpret 
[2] as a single vector space W, then the usual formalisms of bras and kets still 
pass over from the cobordism category. 



^{|^l><0l}(8) 1 = p 



;^CZi^ 10{|^^><0|} = Q 

-; — ^ ^{l0<0|}{lf2>(8)l} = R 



c 



c 



::> = • • R = 1 



PQP = P 

Figure 14 - The Basic Temperley-Lieb Relation 



38 



0> 





o 




Figure 15 - The Key to Teleportation 

Figure 15 illustrates the staightcning of |0) and and the straightening 
of a composition of these applied to \ip), resulting in In the left-hand 
part of the bottom of Figure 15 we illustrate the preparation of the tensor 
product |0) ® \4>) followed by a successful measurement by {n\ in the second 
two tensor factors. The resulting single qubit state, as seen by straightening, 
is \(f)) = OoQlV'). 

From this, we see that it is possible to reversibly, indeed unitarily, transform 
a state via a combination of preparation and measurement just so long as 
the straightenings of the preparation and measurement (0 and Q) are each 
invertible (unitary). This is the key to teleportation [48, 17, 18]. In the 
standard teleportation procedure one chooses the preparation to be (up to 
normalization) the 2 dimensional identity matrix so that \6) = |00) + |11). If the 
successful measurement ft is also the identity, then the transmitted state |0) 
will be equal to j^')- In general we will have ](/)) — . One can then choose a 



39 



basis of measurements each corresponding to a unitary transformation Q 
so that the recipient of the transmission can rotate the result by the inverse of 
ft to reconsitute if he is given the requisite information. This is the basic 
design of the teleportation procedure. 

There is much more to say about the category Co6[0] and its relationship 
with quantum mechanics. We will stop here, and invite the reader to explore 
further. Later in this paper, we shall use these ideas in formulating our rep- 
resentations of the braid group. For now, we point out how things look as we 
move upward to Co6[n] for n > 0. In Figure 16 we show typical cobordisms 
(morphisms) in Co6[l] from two circles to one circle and from one circle to two 
circles. These are often called "pairs of pants" . Their composition is a surface 
of genus one seen as a morphism from two circles to two circles. The bottom 
of the figure indicates a ket-bra in this dimension in the form of a mapping 
from one circle to one circle as a composition of a cobordism of a circle to the 
empty set and a cobordism from the empty set to a circle (circles bounding 
disks). As we go to higher dimensions the structure of cobordisms becomes 
more interesting and more complicated. It is remarkable that there is so much 
structure in the lowest dimensions of these categories. 




Figure 16 - Corbordisms of 1-Manifolds eire Surfaces 
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6 Braiding and Topological Quantum Field The- 
ory 

The purpose of this section is to discuss in a very general way how braid- 
ing is related to topological quantum field theory. In the section to follow, 
we will use the Temperley-Lieb recoupling theory to produce specfic unitary 
representations of the Artin braid group. 

The ideas in the subject of topological quantum field theory (TQFT) are 
well expressed in the book [5] by Michael Atiyah and the paper [78] by Edward 
Witten. Here is Atiyah's definition: 

Definition. A TQFT in dimension d is a functor ^(S) from the cobordism 
category Co6[d] to the category Vect of vector spaces and linear mappings 
which assigns 

1. a finite dimensional vector space Z^E) to each compact, oriented d- 
dimensional manifold E, 

2. a vector Z{Y) G Z{T,) for each compact, oriented {d + l)-dimensional 
manifold Y with boundary E. 

3. a linear mapping Z{Y) : ^(Si) — > Z{Tj2) when Y is a {d+ l)-manifold 
that is a cobordism between Si and S2 (whence the boundary of Y is 
the union of Ei and — E2. 

The functor satisfies the following axioms. 

1. Z{T,^) — Z{I!,y where E^ denotes the manifold E with the opposite 
orientation and ^(E)^ is the dual vector space. 

2. Z(Ei U E2) = ^(Ei) ® ^(E2) where U denotes disjoint union. 

3. If Yi is a cobordism from Ei to E2, I2 is a cobordism from E2 to E3 and 
Y is the composite cobordism Y = Yi Usj Y2, then 

Z(Y) = Z(Y2) o Z(Y,) : Z(Ei) Z(E2) 

is the composite of the corresponding linear mappings. 

4. Z{(j)) — C (C denotes the complex numbers) for the empty manifold (j). 

5. With E X 7 (where I denotes the unit interval) denoting the identity 
cobordism from E to E, Z(E x I) is the identity mapping on Z{Y). 
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Note that, in this view a TQFT is basically a functor from the cobordism 
categories defined in the last section to Vector Spaces over the complex num- 
bers. We have already seen that in the lowest dimensional case of cobordisms 
of zero-dimensional manifolds, this gives rise to a rich structure related to 
quatum mechanics and quantum information theory. The remarkable fact is 
that the case of three-dimensions is also related to quantum theory, and to 
the lower- dimensional versions of the TQFT. This gives a significant way to 
think about three-manifold invariants in terms of lower dimensional patterns 
of interaction. Here follows a brief description. 

Regard the three-manifold as a union of two handlebodies with boundary 
an orientable surface Sg of genus g. The surface is divided up into trinions as 
illustrated in Figure 17. A trinion is a surface with boundary that is topo- 
logically equivalent to a sphere with three punctures. The trinion constitutes, 
in itself a cobordism in Co6[l] from two circles to a single circle, or from a 
single circle to two circles, or from three circles to the empty set. The pattern 
of a trinion is a trivalent graphical vertex, as illustrated in Figure 17. In that 
figure we show the trivalent vertex graphical pattern drawn on the surface of 
the trinion, forming a graphical pattern for this combordism. It should be 
clear from this figure that any cobordism in Co6[l] can be diagrammed by a 
trivalent graph, so that the category of trivalent graphs (as morphisms from 
ordered sets of points to ordered sets of points) has an image in the cate- 
gory of cobordisms of compact one-dimensional manifolds. Given a surface S 
(possibly with boundary) and a decomposition of that surface into triions, we 
associate to it a trivalent graph G{S, t) where t denotes the particular trinion 
decomposition. 

In this correspondence, distinct graphs can correspond to topologically 
identical cobordisms of circles, as illustrated in Figure 19. It turns out that 
the graphical structure is important, and that it is extraordinarily useful to 
articulate transformations between the graphs that correspond to the home- 
omorphisms of the corresponding surfaces. The beginning of this structure is 
indicated in the bottom part of Figure 19. 

In Figure 20 we illustrate another feature of the relationship betweem sur- 
faces and graphs. At the top of the figure we indicate a homeomorphism 
between a twisted trinion and a standard trinion. The homeomorphism leaves 
the ends of the trinion (denoted A,B and C) fixed while undoing the internal 
twist. This can be accomplished as an ambient isotopy of the embeddings in 
three dimensional space that are indicated by this figure. Below this isotopy 
we indicate the corresponding graphs. In the graph category there will have 
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to be a transformation between a braided and an unbraided trivalent vertex 
that corresponds to tliis homeomorphism. 




Figure 17 - Decomposition of a Surface into Trinions 




c 
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Figure 18 - Trivalent Vectors 
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Figure 20 - Tube Twist 



Prom the point of view that we shall take in this paper, the key to the 
mathematical structure of three-dimensional TQFT lies in the trivalent graphs, 
including the braiding of grapical arcs. We can think of these braided graphs 
as representing idealized Fcynman diagrams, with the trivalent vertex as the 
basic particle interaction vertex, and the braiding of lines representing an in- 
teraction resulting from an exchange of particles. In this view one thinks of 
the particles as moving in a two-dimensional medium, and the diagrams of 
braiding and trivalent vertex interactions as indications of the temporal events 
in the system, with time indicated in the direction of the morphisms in the 
category. Adding such graphs to the category of knots and links is an exten- 
sion of the tangle category where one has already extended braids to allow any 
embedding of strands and circles that start in n ordered points and end in 
m ordered points. The tangle category includes the braid category and the 
Temperley-Licb category. These are both included in the category of braided 
trivalent graphs. 

Thinking of the basic trivalent vertex as the form of a particle interaction 
there will be a set of particle states that can label each arc incident to the 
vertex. In Figure 18 we illustrate the labeling of the trivalent graphs by such 
particle states. In the next two sections we will see specific rules for labeling 
such states. Here it suffices to note that there will be some restrictions on these 
labels, so that a trivalent vertex has a set of possible labehngs. Similarly, any 
trivalent graph will have a set of admissible labelings. These are the possible 
particle processes that this graph can support. We take the set of admissible 
labelings of a given graph G as a basis for a vector space V{G) over the 
complex numbers. This vector space is the space of processes associated with 
the graph G. Given a surface S and a decomposition t of the surface into 
trinions, wc have the associated graph G{S, t) and hence a vector space of 
processes V{G{S,t)). It is desirable to have this vector space independent of 
the particular decomposition into trinions. If this can be accomplished, then 
the set of vector spaces and linear mappings associated to the surfaces can 
consitute a functor from the category of cobordisms of one-manifolds to vector 
spaces, and hence gives rise to a one-dimensional topological quantum field 
theory. To this end we need some properties of the particle interactions that 
will be described below. 

A spin network is, by definition a lableled trivalent graph in a category of 
graphs that satisfy the properties outlined in the previous paragraph. We shall 
detail the requirements below. 
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The simplest case of this idea is C. N. Yang's original interpretation of the 
Yang-Baxter Equation [80]. Yang articulated a quantum field theory in one 
dimension of space and one dimension of time in which the i?-matrix giving the 
scattering ampitudes for an interaction of two particles whose (let us say) spins 
corresponded to the matrix indices so that is the amplitude for particles of 
spin a and spin b to interact and produce particles of spin c and d. Since these 
interactions are between particles in a line, one takes the convention that the 
particle with spin a is to the left of the particle with spin b, and the particle with 
spin c is to the left of the particle with spin d. If one follows the concatenation 
of such interactions, then there is an underlying permutation that is obtained 
by following strands from the bottom to the top of the diagram (thinking of 
time as moving up the page). Yang designed the Yang-Baxter equation for R 
so that the amplitudes for a composite process depend only on the underlying 
permutation corresponding to the process and not on the individual sequences 
of interactions. 

In taking over the Yang-Baxter equation for topological purposes, we can 
use the same interpretation, but think of the diagrams with their under- and 
over-crossings as modeling events in a spacetime with two dimensions of space 
and one dimension of time. The extra spatial dimension is taken in displacing 
the woven strands perpendicular to the page, and allows us to use braiding 
operators R and R~^ as scattering matrices. Taking this picture to heart, one 
can add other particle properties to the idealized theory. In particular one 
can add fusion and creation vertices where in fusion two particles interact to 
become a single particle and in creation one particle changes (decays) into two 
particles. These are the trivalent vertices discussed above. Matrix elements 
corresponding to trivalent vertices can represent these interactions. See Figure 
21. 




Figure 21 -Creation and Fusion 

Once one introduces trivalent vertices for fusion and creation, there is the 
question how these interactions will behave in respect to the braiding operators. 
There will be a matrix expression for the compositions of braiding and fusion or 
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creation as indicated in Figure 22. Here we will restrict ourselves to showing the 
diagrammatics with the intent of giving the reader a flavor of these structures. 
It is natural to assume that braiding intertwines with creation as shown in 
Figure 24 (similarly with fusion) . This intertwining identity is clearly the sort 
of thing that a topologist will love, since it indicates that the diagrams can 
be interpreted as embeddings of graphs in three-dimensional space, and it fits 
with our interpretation of the vertices in terms of trinions. Figure 22 illustrates 
the Yang-Baxter equation. The intertwining identity is an assumption like the 
Yang-Baxter equation itself, that simplifies the mathematical structure of the 
model. 





R(g)l 



l(2)R 



R(g)l 
l(g)R 




l(g)R 
R(g)l 

I® R 



Figure 22 - YangBaxterEquation 




Figure 23 - Braiding 
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Figure 24 - Intertwining 



It is to be expected that there will be an operator that expresses the re- 
coupling of vertex interactions as shown in Figure 25 and labeled by Q. This 
corresponds to the associativity at the level of trinion combinations shown in 
Figure 19. The actual formahsm of such an operator will parallel the mathe- 
matics of recoupling for angular momentum. See for example [36]. If one just 
considers the abstract structure of recoupling then one sees that for trees with 
four branches (each with a single root) there is a cycle of length five as shown 
in Figure 26. One can start with any pattern of three vertex interactions and 
go through a sequence of five recouplings that bring one back to the same 
tree from which one started. It is a natural simplifying axiom to assume that 
this composition is the identity mapping. This axiom is called the pentagon 
identity. 




Figure 25 - Recoupling 
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Figure 26 - Pentagon Identity 



Finally there is a hexagonal cycle of interactions between braiding, recou- 
pling and the intertwining identity as shown in Figure 27. One says that the 
interactions satisfy the hexagon identity if this composition is the identity 




Figure 27 - Hexagon Identity 

A graphical three-dimensional topological quantum field theory is an algebra 
of interactions that satisfies the Yang-Baxter equation, the intertwining iden- 
tity, the pentagon identity and the hexagon identity There is not room in this 
summary to detail the way that these properties fit into the topology of knots 
and three-dimensional manifolds, but a sketch is in order. For the case of topo- 
logical quantum field theory related to the group SU (2) there is a construction 
based entirely on the combinatorial topology of the bracket polynomial (See 
Sections 7,9 and 10 of this article.). See [41, 36] for more information on this 
approach. 
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Now return to Figure 17 where we illustrate trinions, shown in relation 
to a trivalent vertex, and a surface of genus three that is decomposed into 
four trinions. It turns out that the vector space V{Sg) ~ V{G{Sg,t)) to 
a surface with a trinion decomposition as t described above, and defined in 
terms of the graphical topological quantum field theory, does not depend upon 
the choice of trinion decomposition. This independence is guaranteed by the 
braiding, hexagon and pentagon identities. One can then associate a well- 
defined vector \M) in V{Sg) whenenver M is a three manifold whose boundary 
is Sg. Furthermore, if a closed three- manifold is decomposed along a surface 
Sg into the union of M_ and M+ where these parts are otherwise disjoint three- 
manifolds with boundary Sg, then the inner product /(M) = (M_|M+) is, up 
to normalization, an invariant of the three- manifold M3. With the definition 
of graphical topological quantum field theory given above, knots and links can 
be incorporated as well, so that one obtains a source of invariants I{M^, K) 
of knots and links in orientable three-manifolds. Here we see the uses of the 
relationships that occur in the higher dimensional cobordism categories, as 
descirbed in the previous section. 

The invariant 7(M^, K) can be formally compared with the Witten [78] integral 

It can be shown that up to limits of the heuristics, Z{M, K) and I{M^, K) are 
essentially equivalent for appropriate choice of gauge group and corresponding 
spin networks. 

By these graphical reformulations, a three-dimensional TQFT is, at base, 
a highly simplified theory of point particle interactions in 2 -|- 1 dimensional 
spacetime. It can be used to articulate invariants of knots and links and 
invariants of three manifolds. The reader interested in the SU{2) case of this 
structure and its implications for invariants of knots and three manifolds can 
consult [36, 41, 59, 19, 64]. One expects that physical situations involving 2-1-1 
spacetime will be approximated by such an idealized theory. There are also 
applications to 3 -|- 1 quantum gravity [6, 7, 50]. Aspects of the quantum Hall 
effect may be related to topological quantum field theory [77]. One can study 
a physics in two dimensional space where the braiding of particles or collective 
excitations leads to non-trival representations of the Artin braid group. Such 
particles are called Anyons. Such TQFT models would describe applicable 
physics. One can think about applications of anyons to quantum computing 
along the fines of the topoological models described here. 
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Figure 28 - A More Complex Braiding Operator 

A key point in the application of TQFT to quantum information theory 
is contained in the structure illustrated in Figure 28. There we show a more 
complex braiding operator, based on the composition of recoupling with the 
elementary braiding at a vertex. (This structure is implicit in the Hexagon 
identity of Figure 27.) The new braiding operator is a source of unitary rep- 
resentations of braid group in situations (which exist mathematically) where 
the recoupling transformations are themselves unitary. This kind of pattern is 
utilized in the work of Freedman and collaborators [24, 25, 26, 27, 28] and in 
the case of classical angular momentum formalism has been dubbed a "spin- 
network quantum simlator" by Rasetti and collaborators [61, 62]. In the next 
section we show how certain natural deformations [36] of Penrose spin net- 
works [66] can be used to produce these unitary representations of the Artin 
braid group and the corresponding models for anyonic topological quantum 
computation. 



7 Spin Networks and Temperley-Lieb Recou- 
pling Theory 

In this section we discuss a combinatorial construction for spin networks that 
generalizes the original construction of Roger Penrose. The result of this gen- 
eralization is a structure that satisfies all the properties of a graphical TQFT 
as described in the previous section, and specializes to classical angular mo- 
mentum recoupling theory in the limit of its basic variable. The construction 
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is based on the properties of the bracket polynomial (as already described in 
Section 4). A complete description of this theory can be found in the book 
"Temperley-Lieb Recouphng Theory and Invariants of Three- Manifolds" by 
Kauffman and Lins [36]. 



The "g-deformed" spin networks that we construct here are based on the 
bracket polynomial relation. View Figure 29 and Figure 30. 
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Figure 29 - Basic Projectors 
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Figure 30 - Two Strand Projector 
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Figure 31 -Vertex 

In Figure 29 we indicate how the basic projector (symmetrizer, Jones- Wenzl 
projector) 
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is constructed on the basis of the bracket polynomial expansion. In this tech- 
nology a symmetrizer is a sum of tangles on n strands (for a chosen integer n) . 
The tangles are made by summing over braid lifts of permutations in the sym- 
metric group on n letters, as indicated in Figure 29. Each elementary braid is 
then expanded by the bracket polynomial relation as indicated in Figure 29 so 
that the resulting sum consists of flat tangles without any crossings (these can 
be viewed as elements in the Tempcrlcy-Lieb algebra). The projectors have the 
property that the concatenation of a projector with itself is just that projector, 
and if you tie two lines on the top or the bottom of a projector together, then 
the evaluation is zero. This general definition of projectors is very useful for 
this theory. The two-strand projector is shown in Figure 30. Here the formula 
for that projector is particularly simple. It is the sum of two parallel arcs and 
two turn-around arcs (with coefficient — l/d, with d — —A^ — A~'^ is the loop 
value for the bracket polynomial. Figure 30 also shows the recursion formula 
for the general projector. This recursion formula is due to Jones and Wenzl 
and the projector in this form, developed as a sum in the Temperley-Lieb 
algebra (see Section 5 of this paper), is usually known as the Jones-Wenzl 
projector. 

The projectors are combinatorial analogs of irreducible representations of a 
group (the original spin nets were based on SU (2) and these deformed nets are 
based on the corresponding quantum group to SU(2)). As such the reader can 
think of them as "particles" . The interactions of these particles are governed 
by how they can be tied together into three- vertices. See Figure 31. In Figure 
31 we show how to tic three projectors, of a, h, c strands respectively, together 
to form a three-vertex. In order to accomplish this interaction, we must share 
lines between them as shown in that Figure so that there are non-negative 
integers i,j,k so that a = i + j,b = j + k,c — i + k. This is equivalent to the 
condition that a + 6 + c is even and that the sum of any two of a, b, c is greater 
than or equal to the third. For example a + b > c. One can think of the vertex 
as a possible particle interaction where [a] and [b] interact to produce [c] . That 
is, any two of the legs of the vertex can be regarded as interacting to produce 
the third leg. 

There is a basic orthogonality of three vertices as shown in Figure 32. Here 
if we tie two three- vertices together so that they form a "bubble" in the middle, 
then the resulting network with labels a and b on its free ends is a multiple of 
an a-line (meaning a line with an a-projector on it) or zero (if a is not equal to 
b). The multiple is compatible with the results of closing the diagram in the 
equation of Figure 32 so the two free ends are identified with one another. On 
closure, as shown in the Figure, the left hand side of the equation becomes a 
Theta graph and the right hand side becomes a multiple of a "delta" where 
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denotes the bracket polynomial evaluation of the a-strand loop with a projector 
on it. The G(a, b, c) denotes the bracket evaluation of a theta graph made from 
three trivalent vertices and labeled with a, b, c on its edges. 

There is a recoupling formula in this theory in the form shown in Figure 33. 
Here there are "6-j symbols" , recoupling coefficients that can be expressed, as 
shown in Figure 35, in terms of tetrahedral graph evaluations and theta graph 
evaluations. The tetrahedral graph is shown in Figure 34. One derives the 
formulas for these coefficients directly from the orthogonality relations for the 
trivalent vertices by closing the left hand side of the recoupling formula and 
using orthogonality to evaluate the right hand side. This is illustrated in Figure 
35. 
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Figure 35 - Tetrahedron Formula for Recoupling Coefficients 
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Finally, there is the braiding relation, as illustrated in Figure 36. 




X' = x(x+2) 

Figure 36 - Local Braiding Formula 

With the braiding relation in place, this g-deformed spin network theory 
satisfies the pentagon, hexagon and braiding naturality identities needed for 
a topological quantum field theory. All these identities follow naturally from 
the basic underlying topological construction of the bracket polynomial. One 
can apply the theory to many different situations. 



7.1 Evaluations 

In this section we discuss the structure of the evaluations for A„ and the 
theta and tetrahedral networks. We refer to [36] for the details behind these 
formulas. Recall that A„ is the bracket evaluation of the closure of the n- 
strand projector, as illustrated in Figure 32. For the bracket variable A, one 
finds that 

One sometimes writes the quantum integer 

A2n 

[n] = i-ir-'A^^. = • 

If 

A = e'""/^^ 
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where r is a positive integer, then 

smyK/r) 

Here the corresponding quantum integer is 

sin{m:/r) 



\n\ 



sin{Tr/r) 



Note that [n + 1] is a positive real number for n = 0, 1,2, ...r — 2 and that 
[r - 1] = 0. 

The evaluation of the theta net is expressed in terms of quantum integers 
by the formula 

R(n h A = ( nm+n+p h + ^ + P+ll'WH'b]! 

^ ' ' ^ ^ ^ [m + n]\[n + p]\\p + m]\ 

where 

a = m + p,b = m + n,c = n + p. 

Note that 

{a + b + c) /2 — m + n + p. 

When A = e*'^/^'', the recoupling theory becomes finite with the restriction 
that only three-vertices (labeled with a, b, c) are admissible when a + b + c < 
2r — 4. All the summations in the formulas for recoupling are restricted to 
admissible triples of this form. 



7.2 Symmetry and Unitarity 

The formula for the recoupling coefficients given in Figure 35 has less symmetry 
than is actually inherent in the structure of the situation. By multiplying all 
the vertices by an appropriate factor, we can reconfigure the formulas in this 
theory so that the revised recoupling transformation is orthogonal, in the sense 
that its transpose is equal to its inverse. This is a very useful fact. It means 
that when the resulting matrices are real, then the recoupling transformations 
are unitary. We shall see particular applications of this viewpoint later in the 
paper. 
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Figure 37 illustrates this modification of the three-vertex. Let Vert[a, b, c] 
denote the original 3-vertex of the Temperley-Lieb recoupling theory. Let 
ModVert[a, h, c] denote the modified vertex. Then we have the formula 



ModVert[a, b, c] = v ^ ° Vert[a, b, c]. 
^jQ{a,b, c) 

Lemma. For the bracket evaluation at the root of unity A = e*'^/^'' the factor 



f{a,b,c) = 



0(a, 6, c) 



is real, and can be taken to be a positive real number for {a, b, c) admissible 
(i.e. a + b + c <2r - 4). 

Proof. By the results from the previous subsection, 

e(a, b,c) = (-l)('^+''+'=)/2e(a, b,c) 

where Q{a,b, c) is positive real, and 

A„A,Ae = (-l)('^+''+<=)[a+ l][c+ 1] 

where the quantum integers in this formula can be taken to be positive real. 
It follows from this that 



/(a, b, c) = 



\ e(a,6,c) 

showing that this factor can be taken to be positive real. □ 

In Figure 38 we show how this modification of the vertex affects the non- 
zero term of the orthogonality of trivalent vertices (compare with Figure 32). 
We refer to this as the "modified bubble identity." The coefficient in the mod- 
ified bubble identity is 



(_l)(6+c-o)/2^ 



[6 + l][c- 



\ [a+1] 

where (a, b, c) form an admissible triple. In particular b + c — a is even and 
hence this factor can be taken to be real. 
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We rewrite the recoupling formula in this new basis and emphasize that 
the recouphng coefficients can be seen (for fixed external labels a, b, c, d) as a 
matrix transforming the horizontal "double- y" basis to a vertically disposed 
double-y basis. In Figures 39, 40 and 41 we have shown the form of this 
transformation,using the matrix notation 

M[a, 6, c, d]ij 

for the modified recoupling coefficients. In Figure 39 we derive an explicit 
formula for these matrix elements. The proof of this formula follows directly 
from trivalent-vertex orthogonality (See Figures 32 and 35.), and is given in 
Figure 39. The result shown in Figure 39 and Figure 40 is the following formula 
for the recoupling matrix elements. 



M[a, 6, c, d]ij = ModTet 



a b i 
c d j 



where vAoAbAcAd is short-hand for the product 



^ ( _ 1 ) (a+6-j)/2 ( _ X ) (c+d-j)/2 ^ _ ^ 



a + l][b+l] 



c+l][d + l]^^^^ 



\ b + \ b + i] 

= {-lY''+'+''+'^/y[a + l][b + l][c+l][d+l] 

In this form, since (a, b,j) and (c, d,j) are admissible triples, we see that this 
coeffient can be taken to be real, and its value is independent of the choice of 
i and j. The matrix M[a, b, c, d] is real-valued. 

It follows from Figure 33 (turn the diagrams by ninety degrees) that 

M[a,b,c,d]~^ = M[b,d,a,c]. 
In Figure 42 we illustrate the formula 

M[a,b, c, df = M[b,d,a,c]. 
It follows from this formula that 

M[a, b, c, d\'^ = M[a, b, c, d\~^. 
Hence M[a,b, c, d] is an orthogonal, real-valued matrix. 
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AaA^Ac 



0( a , b , c ) 



Figure 37 - Modified Three Vertex 
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Figure 39 - Derivation of Modified Recoupling Coefficients 
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Figure 40 - Modified Recoupling Formula 




M[a,b,c,d].. = 

Figure 



a b 
c d 



l-^(iiified Recoupling Matrix 



63 



AaAbAcAd 



a b 


T 


a b 


c d 




c d 



Figure 42 - Modified Matrix Transpose 



Theorem 2. In the Temperley-Lieb theory we obtain unitary (in fact real 
orthogonal) recoupling transformations when the bracket variable A has the 
form A — e*'^/^'' for r a positive integer. Thus we obtain families of unitary 
representations of the Artin braid group from the recoupling theory at these 
roots of unity. 

Proof. The proof is given the discussion above. □ 

In Section 9 we shall show explictly how these methods work in the case of 
the Fibonacci model where A — e^*'^/^. 



8 Fibonacci Particles 

In this section and the next we detail how the Fibonacci model for anyonic 
quantum computing [56, 67] can be constructed by using a version of the two- 
stranded bracket polynomial and a generalization of Penrose spin networks. 
This is a fragment of the Temperly-Lieb recoupling theory [36]. We already 
gave in the preceding sections a general discussion of the theory of spin net- 
works and their relationship with quantum computing. 
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The Fibonacci model is a TQFT that is based on a single "particle" with 
two states that we shall call the marked state and the unmarked state. The 
particle in the marked state can interact with itself either to produce a single 
particle in the marked state, or to produce a single particle in the unmarked 
state. The particle in the unmarked state has no influence in interactions (an 
unmarked state interacting with any state 5* yields that state S). One way 
to indicate these two interactions symbolically is to use a box,for the marked 
state and a blank space for the unmarked state. Then one has two modes of 
interaction of a box with itself: 

1. Adjacency: |^ |^ 
and 



2. Nesting: 



With this convention we take the adjacency interaction to yield a single box, 
and the nesting interaction to produce nothing: 



We take the notational opportunity to denote nothing by an asterisk (*). The 
syntatical rules for operating the asterisk are Thus the asterisk is a stand-in 
for no mark at all and it can be erased or placed wherever it is convenient to 
do so. Thus 



P P p p 



Figure 43 - Fibonacci Particle Interaction 

We shall make a recoupling theory based on this particle, but it is worth 
noting some of its purely combinatorial properties first. The arithmetic of 
combining boxes (standing for acts of distinction) according to these rules 
has been studied and formalized in [73] and correlated with Boolean algebra 
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and classical logic. Here within and next to are ways to refer to the two 
sides delineated by the given distinction. From this point of view, there are 
two modes of relationship (adjacency and nesting) that arise at once in the 
presence of a distinction. 



P P 




dim(vj'''') 



= 1 









Figure 44 - Fibonacci Trees 



Prom here on we shall denote the Fibonacii particle by the letter P. Thus 
the two possible interactions of P with itself are as follows. 

1. P,P — >* 

2. P,P — >P 

In Figure 43 we indicate in small tree diagrams the two possible interactions 
of the particle P with itself. In the first interaction the particle vanishes, 
producing the asterix. In the second interaction the particle a single copy of 
P is produced. These are the two basic actions of a single distinction relative 
to itself, and they constitute our formalism for this very elementary particle. 

In Figure 44, we have indicated the different results of particle processes 
where we begin with a left-associated tree structure with three branches, all 
marked and then four branches all marked. In each case we demand that the 
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particles interact successively to produce an unmarked particle in the end, at 
the root of the tree. More generally one can consider a left-associated tree 
with n upward branches and one root. Let T{ai, 02, • • • , a„ : 6) denote such a 
tree with particle labels Oi, • • • , a„ on the top and root label b at the bottom of 
the tree. We consider all possible processes (sequences of particle interactions) 
that start with the labels at the top of the tree, and end with the labels at 
the bottom of the tree. Each such sequence is regarded as a basis vector in a 
complex vector space 

associated with the tree. In the case where all the labels are marked at the 
top and the bottom label is unmarked, we shall denote this tree by 

T/lll-ll _ T/(n) 
Vq — Vq 

where n denotes the number of upward branches in the tree. We see from 
Figure 44 that the dimension of Vq is 1, and that 

dim{vi^^) = 2. 

This means that Vo^^^ is a natural candidate in this context for the two-qubit 
space. 

Given the tree r(l, 1, 1, ■ ■ • , 1 : 0) (n marked states at the top, an unmarked 
state at the bottom), a process basis vector in V^*^"^ is in direct correspondence 
with a string of boxes and asterisks (I's and O's) of length n—2 with no repeated 
asterisks and ending in a marked state. See Figure 44 for an illustration of the 
simplest cases. It follows from this that 

where fk denotes the k-th Fibonacci number: 

/0 = l,/i = l,/2 = 2,/3 = 3,/4 = 5,/5 = 8,--- 

where 

fn+2 — fn+1 + fn- 

The dimension formula for these spaces follows from the fact that there are fn 
sequences of length n — 1 of marked and unmarked states with no repetition 
of an unmarked state. This fact is illustrated in Figure 45. 
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PPP PR* P*P *pp*p^ 
Tree of squences with no occurence of * * 

Figure 45 - Fibonacci Sequence 



9 The Fibonacci Recoupling Model 



We now show how to make a model for recouphng the Fibonacci particle by 
using the Temperley Lieb recoupling theory and the bracket polynomial. Ev- 
erything we do in this section will be based on the 2-projector, its properties 
and evaluations based on the bracket polynomial model for the Jones poly- 
nomial. While we have outlined the general recoupling theory based on the 
bracket polynomial in earlier sections of this paper, the present section is self- 
contained, using only basic information about the bracket polyonmial, and the 
essential properties of the 2-projcctor as shown in Figure 46. In this figure we 
state the definition of the 2-projector, list its two main properties (the opera- 
tor is idempotent and a self-attached strand yields a zero evaluation) and give 
diagrammatic proofs of these properties. 
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a = 



-1/5 



u 

n 



ci = 



u 



-1/5 



u_ 





-(1/5)5 



= 



-1/8 



n 



Figure 46 - The 2-Projector 



In Figure 47, we sliow tlie essence of tlie Temperley-Lieb recoupling model 
for the Fibonacci particle. The Fibonaccie particle is, in this mathematical 
model, identified with the 2-projector itself. As the reader can see from Figure 
47, there are two basic interactions of the 2-projector with itself, one giving 
a 2-projector, the other giving nothing. This is the pattern of self-iteraction 
of the Fibonacci particle. There is a third possibility, depicted in Figure 47, 
where two 2-projcctors interact to produce a 4-projector. We could remark at 
the outset, that the 4-projector will be zero if we choose the bracket polynomial 
variable A — e''^!^ . Rather than start there, we will assume that the 4-projector 
is forbidden and deduce (below) that the theory has to be at this root of unity. 
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Figure 47 - Fibonacci Pcirticle as 2-Projector 

Note that in Figure 47 we have adopted a single strand notation for the particle 
interactions, with a solid strand corresponding to the marked particle, a dotted 
strand (or nothing) corresponding to the unmarked particle. A dark vertex 
indicates either an interaction point, or it may be used to indicate the single 
strand is shorthand for two ordinary strands. Remember that these are all 
shorthand expressions for underlying bracket polynomial calculations. 

In Figures 48, 49, 50, 51, 52 and 53 we have provided complete diagram- 
matic calculations of all of the relevant small nets and evaluations that are 
useful in the two-strand theory that is being used here. The reader may wish 
to skip directly to Figure 54 where we determine the form of the recoupling 
coefficients for this theory. We will discuss the resulting algebra below. 

For the reader who does not want to skip the next collection of Figures, 
here is a guided tour. Figure 48 illustrates three three basic nets in case of 
two strands. These are the theta, delta and tetrahedron nets. In this Figure 
we have shown the decomposition on the theta and delta nets in terms of 2- 
projectors. The Tetrahedron net will be similarly decomposed in Figures 52 
and 53. The theta net is denoted ©, the delta by A, and the tetrahedron by T. 
In Figure 49 we illustrate how a pedant loop has a zero evaluation. In Figure 
50 we use the identity in Figure 49 to show how an interior loop (formed by 
two trivalent vertices) can be removed and replaced by a factor of G/A. Note 
how, in this figure, line two proves that one network is a multiple of the other, 
while line three determines the value of the multiple by closing both nets. 
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Figure 51 illustrates the explicit calculation of the delta and theta nets. The 
figure begins with a calculation of the result of closing a single strand of the 
2- projector. The result is a single stand multiplied by {5 — 1/5) where 5 — 
—A^ — A~'^, and A is the bracket polynomial parameter. We then find that 

A = 5^ - 1 

and 

e = (5 - l/5f5 - A/5 ^{5- l/5){5'' - 2). 

Figures 52 and 53 illustrate the calculation of the value of the tetrahedral 
network T. The reader should note the first line of Figure 52 where the tetrad- 
edral net is translated into a pattern of 2-projectors, and simplified. The rest 
of these two figures are a diagrammatic calculation, using the expansion for- 
mula for the 2-projector. At the end of Figure 53 we obtain the formula for 
the tetrahedron 

T={5-l/5f{5^-2)-2e/6. 




Figure 48 - Theta, Delta and Tetrahedron 
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Figure 50 - LoopEvaluation-2 
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= (6-1/5)^5 - A/5 



Figure 51 - Calculate Theta, Delta 




Figure 52 - Calculate Tetrahedron — 1 
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Figure 53 - Calculate Tetrahedron - 2 



Figure 54 is the key calculation for this model. In this figure we assume that 
the recoupling formulas involve only and 2 strands, with corresponding to 
the null particle and 2 corresponding to the 2-projector. (2 + 2 = 4 is forbidden 
as in Figure 47.) From this assumption we calculate that the recoupling matrix 
is given by 





Figure 54 - Recoupling for 2-Projectors 
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Figure 55 - Braiding at the Three- Vertex 
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Figure 56 - Braiding at the Null-Three- Vertex 

Figures 55 and 56 work out the exact formulas for the braiding at a three- vertex 
in this theory. When the 3-vertex has three marked lines, then the braiding 
operator is multiplication by — A^, as in Figure 55. When the 3-vertex has two 
marked lines, then the braiding operator is multiplication by A^, as shown in 
Figure 56. 
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Notice that it follows from the symmetry of the diagrammatic recoupling for- 
mulas of Figure 54 that the square of the recoupling matrix F is equal to the 
identity. That is, 

/ 1 A _ 2 _ M/A A/0 W 1/A A/e \ _ 

\^ 1 )~ ~ [ e/A2 TA/e^ ) [ e/A2 ta/o^ ) ~ 

( i/A^ + i/A i/e + TAVe^ \ 
1^ e/A^ + T/(Ae) i/A + A^rve^ ) ■ 

Thus we need the relation 

l/A+l/A^^ 1. 
This is equivalent to saying that 

A^ = 1 + A, 
a quadratic equation whose solutions are 

A = (1± V5)/2. 

Furthermore, we know that 

A = 6^-1 

from Figure 51. Hence 

A^ = A + 1 = 5l 
We shall now specialize to the case where 

A = 5 = (1 + V5)/2, 
leaving the other cases for the exploration of the reader. We then take 

so that 

5 = -A^ - A-"^ = -2cos(67r/5) = (1 + \/5)/2. 
Note that 5-1/5=1. Thus 

e = (5-1/5)^5- A/5 = 5-1. 

and 

T = (5 - 1/5)2(5' - 2) - 20/5 = (5" - 2) - 2(5 - l)/5 
= (5- l)(5-2)/5 = 35-5. 

Note that 

T = -eVA', 

from which it follows immediately that 

This proves that we can satisfy this model when A = 5 = (1 + -\/5)/2. 
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For this specialization we see that the matrix F becomes 

/ 1/A A/e \_( 1/A A/0 \_ / 1/A A/0 \ 

[ 0/A2 TA/02 J [ 0/A2 (-0VA2)A/02 J [ 0/A2 -1/A J 

This version of F has square equal to the identity independent of the value of 
0, so long as A^ = A + 1. 

The Final Adjustment. Our last version of F suffers from a lack of symme- 
try. It is not a symmetric matrix, and hence not unitary. A final adjustment 
of the model gives this desired symmetry. Consider the result of replacing each 
trivalent vertex (with three 2-projector strands) by a multiple by a given quan- 
tity a. Since the has two vertices, it will be multiplied by a^. Similarly, the 
tetradhedron T will be multiplied by «*^. The A and the d will be unchanged. 
Other properties of the model will remain unchanged. The new recoupling 
matrix, after such an adjustment is made, becomes 

( 1/A 

For symmetry we require 
We take 

a' = VA3/0. 

With this choice of a we have 

A/(q;'0) = A0/(0VA3) = 1/VA. 
Hence the new symmetric F is given by the equation 

1/A I/VA 





where A is the golden ratio and r = 1/A. This gives the Fibonacci model. 
Using Figures 55 and 56, we have that the local braiding matrix for the model 
is given by the formula below with A — e^"^'/^. 

_ / \ 

-e''^'/^ ■ 



R 




The simplest example of a braid group representation arising from this 

theory is the representation of the three strand braid group generated by = 
R and 5*2 = FRF (Remember that F = F^ = F^^.). The matrices 5*1 and 5*2 
are both unitary, and they generate a dense subset of the unitary group U{2), 
supplying the first part of the transformations needed for quantum computing. 
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10 Quantum Computation of Colored Jones 
Polynomials and the Witten-Reshetikhin- 
Turaev Invariant 

In this section we make some brief comments on the quantum computation 
of colored Jones polynomials. This material will be expanded in a subsequent 
publication. 




J 



B 





J 






P(B) 








= if b 4 







a, yava-a 





= B(0,0)(Aa) 



Figure 57 - Evaluation of the Plat Closure of a Braid 
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First, consider Figure 57. In that figure we illustrate the calculation of the 
evalutation of the (a) - colored bracket polynomial ioi the plat closure P{B) of 
a braid B. The reader can infer the definition of the plat closure from Figure 
57. One takes a braid on an even number of strands and closes the top strands 
with each other in a row of maxima. Similarly, the bottom strands are closed 
with a row of minima. It is not hard to see that any knot or link can be 
represented as the plat closure of some braid. 




Figure 58 - Dubrovnik Polynomial Specialization at Two Strands 

The (a) - colored bracket polynonmial of a link L, denoted < L >„, is the 
evaluation of that link where each single strand has been replaced by a parallel 
strands and the insertion of Jones- Wenzl projector (as discussed in Section 7). 
We then see that we can use our discussion of the Temperley-Lieb recoupling 
theory as in sections 7,8 and 9 to compute the value of the colored bracket 
polynomial for the plat closure PB. As shown in Figure 57, we regard the 
braid as acting on a process space Vq'"'''"'"' and take the case of the action on 
the vector v whose process space coordinates are all zero. Then the action of 
the braid takes the form 

Bv{Q, • • • , 0) = Ea;^,...,a;„S(xi, • • • , Xn)v{xi, ■■■,Xn) 
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where B{xi, ■ ■ ■ , Xn) denotes the matrix entries for this recouphng transforma- 
tion and v{xi, ■ ■ ■ , Xn) runs over a basis for the space Vq'""''"'"". Here n is even 
and equal to the number of braid strands. In the figure we illustrate with 
n = 4. Then, as the figure shows, when we close the top of the braid action 
to form PB, we cut the sum down to the evaluation of just one term. In the 
general case we will get 

<PB >„=S(0,---,0)A-/2. 

The calculation simplifies to this degree because of the vanishing of loops in 
the recoupling graphs. The vanishing result is stated in Figure 57, and it is 
proved in the case a = 2 in Figure 49. 

The colored Jones polynomials are normalized versions of the colored bracket 
polymomials, differing just by a normalization factor. 

In order to consider quantumn computation of the colored bracket or col- 
ored Jones polynomials, we therefore can consider quantum computation of 
the matrix entries -5(0, ■ ■ ■ , 0). These matrix entries in the case of the roots of 
unity A = e*'^/^'' and for the a = 2 Fibonacci model with A = e^*'^/^ are parts 
of the diagonal entries of the unitary transformation that represents the braid 
group on the process space Vq'°''"''°'. We can obtain these matrix entries by us- 
ing the Hadamard test as described in section 4- As a result we get relatively 
efficient quantum algoritms for the colored Jones polynonmials at these roots 
of unity, in essentially the same framework as we described in section 4, but 
for braids of arbitrary size. The computational complexity of these models is 
essentially the same as the models for the Jones polynomial discussed in [1]. 
We reserve discussion of these issues to a subsequent publication. 

It is worth remarking here that these algorithms give not only quantum 
algorithms for computing the colored bracket and Jones polynomials, but also 
for computing the Witten-Reshetikhin-Turaev (WRT) invariants at the above 
roots of unity. The reason for this is that the WRT invariant, in unnormalized 
form is given as a finite sum of colored bracket polynomials: 

WRT{L) = ^IzlAa < L >a, 

and so the same computation as shown in Figure 57 applies to the WRT. This 
means that we have, in principle, a quantum algorithm for the computation 
of the Witten functional integral [78] via this knot-theoretic combinatorial 
topology. It would be very interesting to understand a more direct approach 
to such a computation via quantum field theory and functional integration. 
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Finally, we note that in the case of the Fibonacci model, the (2)-colored 
bracket polynomial is a special case of the Dubrovnik version of the Kauffman 
polynomial [38]. See Figure 58 for diagammatics that resolve this fact. The 
skein relation for the Dubrovnik polynomial is boxed in this figure. Above the 
box, we show how the double strands with projectors reproduce this relation. 
This observation means that in the Fibonacci model, the natural underlying 
knot polynomial is a special evaluation of the Dubrovnik polynomial, and the 
Fibonacci model can be used to perform quantum computation for the values 
of this invariant. 
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